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ABSTRACT 

Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed 
the Chandrasekhar mass and explode as a type la supernova? We explore the range of accretion rates that 
allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the 
initial mass that will allow it to reach 1 .4Mq — the Chandrasekhar mass. We follow the evolution through 
a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the 
effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. 
We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass 
lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, 
that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during 
subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the 
helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout 
the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady 
state, instead of explosive. We thus And a remarkably large parameter space within which long-term, 
self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite 
its helium flashes. 

Subject headings: binaries: close — novae, cataclysmic variables — supemovae: type la — white dwarfs 


1. Introduction 

The stellar evolutionary path or paths leading to a 
Type la Supernova (SNIa) has been an open question 
for decades. Unveiling the identity of the progeni¬ 
tors of these powerful explosions, which is the stan¬ 
dard candle that demonstrated the existence of dark 
energy, could have important implications for our un¬ 
derstanding of the acceleration of the expansion of the 
universe. The two possible paths that may lead to 
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a SNIa are the double degenerate (DD) and the sin¬ 
gle degenerate (SD) scenarios. The DD scenario is 
the case of merging binary white dwarfs (WD) with a 
combin ed mass that exceeds the Chandrasekha r mass 
(McsA dlben & Tutukovlflosi IWebbinklll984 . The 
SD scenario is the case of a semi-detached binary 
composed of a WD and a secondary, where the WD 
grows in mass by accreting hydrogen at a high rate, 
leading to a long series of rel atively mild nova out 
bursts, until it reaches Mc^n ( Whelan & IbenI 1973 
Hillebrandt & NiemeveilEoOOll . 


The strongest criticisms of the SD scenario have 
been that: (1) hydrogen accretion onto a WD usu- 
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ally l eads to nova erup tions which decrease the WD 
mass ( Yaron et'^l2005 ). and ( 2 ) even when hydrogen 
flashes lead to a growing WD mass, the accumulated 
helium shells must themselves erupt a nd be ejected 
( Idan et al. 2013 : Newsham et al. 2013 ). In this pa¬ 
per we focus on the SD scenario, in an effort to, for 
the first time, self-consistently determine whether long 
term accretion of hydrogen can increase the mass of a 
WD up to the Chandrasekhar limit. In particular, we 
investigate the parameter space within which the WD 
will eject less mass than it has accreted at the end of 
each nova cycle, thus growing toward Mch. 

Thi s paper is a direct co ntinuation of our previous 
work (IHillman et al.ll2015h . in which we have shown 
the range of accretion rates that may “push” a 1 . 4 M 0 
WD toward a SNIa. Here we expand the investigation 
over the entire parameter space of accretion rate (M) 
and the initial WD mass (Mwd) that can lead to Mch- 
Both the Hubble time and the initial secondary mass 
(Mg) provide additional constraints, and we include 
those too. These systems are, by their natures, recur¬ 
rent novae (RN), and we deduce signatures of these 
systems that may be observationally detectable. 

Observatio ns of novae on massiv e WDs have long 
been reported. Hachisu et al. ( 2000l) reproduced a the¬ 
oretical light curve of the 1999 outburst of the recur¬ 
rent nova U Sco and conclude the mass of the WD to be 
1 . 37 ±O.O 1 M 0 . Hachisu & Katol (2002) numerically 
reproduced light curves of the 1998 outburst of V2487 
Oph, and deduced the WD mass to be 1 . 35 ±O.O 1 M 0 
with an accretion rate of ^ 1 . 5 xlO“^M 0 yr“^ and 
a recurrence period of ^40 years. They estimate 
the growth rate to be ^ 2 xlO“®M 0 ?/r“^ and there¬ 
fore conclude that t his system is a strong SNIa pro¬ 
genitor candidate. Zamanov et al.l (120031) estimated 
the mass of the WD in the recurrent nov a system T 
CrB to be ^ 1 . 34 M 0 . Kato et al. (2008) estimated 


the mass of the WD in the recurrent nova system 
RS Oph to be '^ 1 . 35 M 0 and growing at a rate of 
^5x10“® — MQyr~^. iRajoelimanana et al 


(120131) observed a periodicity of ^450 and ~180 days 
in the Cal 83 and RX J0513.9-6951 systems respec¬ 
tively, implying WDs wit h masses of ^1 . 25 — 
and high accretion rates. ISahman et al.l (12013h report 
a WD mass of l.OO±O.14M0 for the RN system Cl 
Aql an d estimate it wi l l evolv e into a SNIa within 10^ 
years. Baneriee et al. ( 20141) analyzed the 2014 out¬ 
burst of the RN system V745 Sco and suggest the WD 
must be highly mas sive and the sy s tem to be a potential 


(120141) and iTang et al.l (120141) have reported observa¬ 
tions of the RN system RX J0045.4H-4154 in the An¬ 
dromeda galaxy, which is estimated to have a WD with 
a mass of at least 1 . 3 M 0 , a recurrence time of year, 
and an accretion rate of at least 1 . 7 xlO“^M 02 /r“^. 

It is widely believed that only very massive WDs, 
accreting at very high rates, can reach Mch- Numeri¬ 
cal simulations have shown that such systems (massive 
WDs and high rates of accretion) eject less mass than 
has be e n accreted during e a ch nov a (IKovetz & Prialnik 


19941: Prialnik & Kovetz 1995: Yaron etal. 2005t 


Starrfleld et al. 2012bl: Hillman et al. 2015 ). thus al¬ 


lowing the WD to grow and rendering it a SNIa pro¬ 
genitor candidate. However, no fully self-consistent 
numerical study of hydrogen and helium accretion, 
over millions of years, has ever been carried out. The 
“belief’ that such systems yield SNIa is just that — a 
belief, not grounded in detailed simulations. The goal 
of this study is to provide the first such extensive and 
self-consistent simulations. 

In the next section, )J2] we describe the parameter 
space determined via hydrogen accretion and recurrent 
nova flashes, followed by our results of long-term he¬ 
lium accretion and flashes in jj3] A discussion and our 
conclusions are presented in S|4] 

2. Hydrogen accretion simnlations 
2.1. Method of calculation 

The calculations were carried out using a hy 


drodyn amic code developed by iPrialnik & Kovetz 


(1995) and described thoroughly in IHillman et al 


(l2015l) . We covered WD masses ranging from 0.65 
to 1 . 4 M 0 and accretion rates ranging from 3x10“® 
to 6 xlO~^M^jvr~ ^. These are the limits on M that 


SNIa progenitor. iDarnlev et al.l (l2014l) . iHenze et al 


Hillman et al.l (120151) have determined to be the range 
in which a WD retains part (or all) of the mass that 
it has accreted during a nova cycle, thus allowing it 
to secularly grow in mass. For each (Mwd.AF) com¬ 
bination, we allowed the code to run over a few tens 
of consecutive nova cycles, in order to understand the 
differences in a WD’s behavior for different parame¬ 
ter combinations. For each simulation, we recorded 
the accreted mass (m-acc) and the ejected mass (mej) 
at the end of each nova cycle. We also recorded the 
duration {D) of each nova cycle, defined as the time 
between successive hydrogen flash peak temperatures 
in the burning shell. We used these data to calculate 
the time (r) needed for the WD to reach Mch. as well 
as the minimal mass of the secondary star required to 
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supply sufficient mass for the WD to reach Mch- 


progressively higher than the current Mwd^ 


2.2. The parameter space of SNIa candidates - 
Results of evolutionary calculations 

For each parameter combination, that is, for each 
{M^Yy,M) pair, we calculated an effective accretion 
rate (Meff). which is the net average accretion rate 
throughout a cycle. 


Meff(MwD,M)= (1) 

FiglUdisplays Mgg, spanning the ranges of Mwd and 
of M. For high values of M, M^e grows with the 
WD mass, and as M is reduced, the trend becomes 
less pronounced. This means that, for the higher ac¬ 
cretion rates, Meff will become even higher as the WD 
grows in mass, thus accelerating the growth rate. The 
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Fig. 1.— The effective accretion rate (Meff) vs. the 
WD mass (Mwd) and the actual hydrogen accretion 
rate (M). 


effective rate of accretion can now be used to calculate 
the time, r, required for a WD of a given mass, accret¬ 
ing hydrogen-rich mass at a given rate, to reach Mch 
and explode as a SNIa: 


r(MwD, M) 


r-Mch 


dm 


'Mwd 


( 2 ) 




M,+i - M, 
-Meff , 


( 3 ) 


where the index i runs over the Mwd series, with 
Mi(,=MwD and Mi^=Mch- The evolution of t is 
shown in Fig|2] where each point expresses the time 
required to reach Mch from the current Mwd- The 
required time is shorter for higher accretion rates, 
not only because the mass is being accreted faster, 
but also because the net accreted mass is larger for 
higher accretion rates, i.e., WDs accreting at higher 
rates lose less of t he ac c reted mass at each cycle 
(Kovetz & Prialnik 19941: Prialnik & Kovetz 19951: 


Yaron et al.l2005l: Starrfield et al.l2012bl: Ftillman et al 


2015h . The reason is that at higher accretion rate s, the 
accretion phase is shorter ( Prialnik et al. 1982h . and 
so is the time allowed for diffusion of hydrogen in¬ 
ward. Therefore, the thermonuclear runaway (TNR) is 
weaker and occurs cl oser to the surface, whi ch results 
in less ejected mass ( Starrfield et alJEoi2ah . Exami¬ 
nation of the top left corner of Fig |2] reveals that the 
least massive WD in the sample (0.65Mq), accret¬ 
ing at rates of 5xl0“®MQ?/r“^ and less, will require 
more than a Hubble time in order to reach Mch- Thus, 
these low-mass, slowly accreting models are not pos¬ 
sible SNIa candidates at the current epoch of cosmic 
history. But remarkably, if the donor star is massive 
enough and M is large enough, even a 0.65Mq WD 
can be grown to the Chandrasekhar mass. 
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Since we do not continuously follow the evolution of 
each model up to the Chandrasekhar mass, the integral 
in eq.dHi is calculated piecewise, using the results of 
runs with the same accretion rate and initial masses 


Fig. 2.— The required time (r) to reach the Chan¬ 
drasekhar mass vs. the current WD mass (Mwd) 
for different rates of accretion (mass transfer) of 
hydrogen-rich matter (M), as shown. 
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We emphasize that in determining whether a WD 
of initial mass Mwd can reach Mch, two important 
constraints are in play. The first, as already noted, a 
WD must be able to grow to Mch within a Hubble 
time. The second is the initial mass of the donor Mg 
which poses a constraint as well. A lower limit for Mg, 
shown in Fig [3 is simply obtained by 

Afs,min(AfwD, AT) = Mt(Mwd, AT) (4) 



-•-3E-8IVI 

Suny^'' 

-.-4E-8IVI 

Suny^'' 

-a^5E-8IVI 

Suny^'' 

— 1E-7IVI 

Suny^'' 

hi,-2E-7IVI 

Suny^'' 

4E-7IVI 

Suny^'' 

-.-5E-7IVI 

Suny^'' 

^6E-7IVI 

Suny^'' 


WD'- Sun-' 


Fig. 3.— Lower limit on the donor mass (ATg) required 
to grow a WD companion to the Chandrasekhar mass 
vs. the WD mass (Mwd), for different rates of ac¬ 
cretion (mass transfer) of hydrogen-rich matter (AT), 
as marked. The black line corresponds to MwD=ATg, 
which defines the upper limit for stable mass transfer, 
assuming Roche-lobe overflow. This limit does not ap¬ 
ply for wind accretion in a symbiotic binary. 


For example, a WD of mass 1.2Mq, accreting at a 
rate of 5xlQ~^M 0 yr~^, will require a donor mass of 
at least ^3Mq to reach Mph, but a donor of ~ 1 Mq 
will suffice, if the accretion rate is about twice as high. 

In the first case, accretion from a red giant wind in 
a symbiotic binary may provide the required condi¬ 
tions, while in the second, accretion in a close binary 
by Roche-lobe overflow is the likely scenario. 

In conclusion, the parameter space where a SNIa 
may still be obtained is bordered by all the masses 
within our range (i.e., O.65 <Mwd< 1-4M0), accret¬ 
ing at rates higher than 5xlO“®MQ?/r“^, and masses 
higher than 1 . 2 M 0 accreting at any rate within the lim¬ 
its we have defined (i.e., 3xlO“®<AT<6x lO“^M 0 yr“^). 


2.3. Implications for observations 

Two important features of the nova cycle are the 
cycle duration, D, and the hydrogen flash duration, /, 
where the latter is a fraction of the former and is de¬ 
fined as the time from the beginning of the luminos¬ 
ity rise until the end of the luminosity decline. For 
the cases considered here, unlike classical novae, / 
is a non-negligible fraction of D, and both are short 
enough to be measurable directly with observations. 
We show D in Figl4]on a double logarithmic scale. 
There is a strong inverse relationship between D and 
both Mwd and M: the cycle becomes drastically 
shorter as the WD mass grows and decreases with in¬ 
creasing accretion rate. 
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Fig. 4.— Cycle duration D (the time between suc¬ 
cessive hydrogen flashes) vs. accretion rate (AT) on 
a double logarithmic scale, for different WD masses, 
with linear fits. 


Hillman et al.l (12015h determined the hydrogen flash 


duration for a range of WD masses, but always ac¬ 
creting at one constant rate of 5 xlO-^M 0 yr-i. We 
have now extended these calculations, and in Fig|3 
we present the flash duration / throughout the entire 
(Mwd,AT) grid. The hydrogen flash duration varies 
immensely with the change in WD mass. The lower 
the WD mass, the longer the flash duration. By con¬ 
trast, the rate of accretion has a minor effect on the 
duration of the flash. Thus, by using the two observ¬ 
ables, D and /, we can derive the key parameters of a 
RN, as follows. The WD mass can be estimated from 
Figl5]based on the observed / alone, via the following 
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Fig. 5.— Hydrogen flash duration (/) vs. WD mass 
for different rates of accretion (mass transfer) 
of hydrogen-rich matter (M), as shown. Note the very 
weak dependence on accretion rate, which results in a 
/(Mwd) relation represented by the dotted line. 



Fig. 6 .— The coefficients, A and B as functions of the 
WD mass (Mwd), corresponding to Eq| 6 ]and Fig|4] 


quadric function fit: 

Mwd = -0.04284(log/)2 -p 0.02884 log / + 1.430 
= 0.9722) 

(5) 

where / is given in days. The estimated Mwd can 
then be used together with the observed D, using the 
results shown in FigH] to estimate the average rate at 
which mass is being accreted. Fig|4]shows a very close 
linear dependence of log D on log M for each value of 
-^WD- This may be inverted to obtain: 


logM = —A(MwD)log71 — i?(MwD) 
(average) = 0.9899) 


where D is in years and the coefficients, A and B, are 
plotted in Figl 6 ]as a function of the WD mass. 


For example, the RN system R S Oph displays a 
flash duration of 115 — 140 days (Katoetal. 2008 : 
Adamakis et al.l Uoi ll: Osborne et al. 2011 ). Based 


on FiglS] and Eql5] this corresponds to a mass of 
^ 1.3 M(7). It is known to have erupted every ^20 
years dZ atnanov et al.l 2010t lAdamakis et al.l l201 it 
Vavtet et al.l201 1 ). According to EiglHand using Eig]^ 


and Eqlb] these values of Mwd and D correspond to 
an accretion rate of ^ SxlQ~^ MQyr~^. These results 
are in agre ement with calcul ations made by other au¬ 
thors, e.g., Kato et al. ( 2008h who have estimated the 
mass of RS Oph to be ^I.SSMq and the growth rate 
to be 5x10“® — lxlO“^M 0 j/'r“^. 


3. Helium accretion simulations 


3.1. Helium flashes in the literature 


The result of recurrent nova cycles, described 
above, is the accumulation of a helium-rich layer on 
top of the initial WD core. In these RN cases, part or 
all of the accreted hydrogen is retained at the end of the 
outburst and is burned into helium. The maximal tem¬ 
perature never rises much above the threshold temper¬ 
ature for helium igniti on for extended periods of time 
(e.g., .lose et all 1993 ). Thus, during these RN erup¬ 
tions, negligible amounts of the accumulating helium 
is fused into carbon or heavier elements. Eventually 
a sufficiently massive helium layer will accumulate, 
and at the bottom of this layer, pressure will build up 
to high enough levels for triggering helium ignition 
and causing a powerful TNR. Many theoretical stud¬ 
ies have reported the occurrence o f powerful helium 
flashes and resultant mass ejection. Ildan et al.l (120 13h 
simulated the accretion of solar composition mate¬ 
rial on WDs with masses in the range 1.0 — 1 . 4 M 0 
accreting at a rate of lO“®M 0 yr“^ and obtained a he¬ 
lium flash after ~4000 hydrogen flashes which ejected 
nearly all of the accreted mass. Eor a 1 . 35 M 0 WD 
accreting solar composition material at rates in the 
range 5X10“^ — 3.2x lO“®M 02 /r“^, Newsham et al. 


(12013h obtained steady hydrogen burning that resulted 
in helium flashes. 

Kato & Hachisu ( 1999li simulated helium accretion 


onto a \ .3Mq WD and obtained helium flashes with a 
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mass accumulation efficiency ranging from 1 to 0.385 
for helium accretion rates ra nging from ^1.26xl0~^ 
to ~ 5 .OxlO“®M 0 j/r“^. Kato & Hachisu ( 2004 1 


simulated helium accretion for WD masses in the 
range 0.7 — 1.35Mq and accretion rates of 10“^ and 
10 ~^MQyr~^, and obtained helium flashes and corre¬ 
sponding accretion efficiency factors implying that in 
all cases a fraction of the accreted helium is retained 


and bu rned into carbon and oxygen. iPiersanti et al 


(120141) simulated helium accretion onto WDs of 
masses in the range 0.6 — 1 .O 2 M 0 and a wide range of 
accretion rates, and obtained helium flashes for accre¬ 
tion rates in the range ^3x10“® —3.5 xlO“®M 0 j/r“^ 
(depending on the WD mass), while the lower rates 
within the regime lead to strong flashes and mass loss 
and the higher rates within the re gime lead t o mild 
flashes and some mass retention. Ijose et alJ (Il993h 


and Cassisi et al.l ( 1998 ) simulated helium flashes both 
via the accretion of solar composition material that 
leads to helium accumulation and via direct helium ac¬ 
cretion. The former studied WDs with low masses of 
O. 5 I 6 M 0 and Q.SMq while the latter considered more 
massive WDs in the range 1.0 — 1.2Mq. Both used 
various accretion rates and obtained helium flashes 
at high accretion rates after a few tens of hydrogen 
flashes. In addition, they both obtained helium flashes 
via direct helium accretion, and reported the time to 
a helium flash to be longer when accreting helium di¬ 
rectly. They found that the amount of mass lost is 
inversely dependent on the accretion rate. 


In these calculations, hydrogen accretion is re¬ 
placed by helium accretion in order to avoid the ex¬ 
tremely time-consuming simulation of thousands or 
tens of thousands of recurrent hydrogen flashes that 
occur between helium flashes. The conclusion of these 
studies is that taking into account the loss of helium- 
rich material in the helium flashes, in addition to the 
hydrogen-rich material ejected in the hydrogen flashes, 
renders the retention of mass negligible and the growth 
rate of the WD extremely slow and inefficient. This 
conclusion invalidates, in principle, the possibility of 
a WD attaining the Chandrasekhar mass by accretion 
from a close companion. However, the above studies 
involving the simulation of helium flashes do not con¬ 
tinue beyond the first few helium flashes, and hence 
their conclusion of negligible mass retention is based 
on an unsubstantiated extrapolation. 


3.2. Long-term evolutionary calculations 

Our next step is therefore to examine the question 
of helium flashes and their long-term consequences. 
We have adopted the common procedure of replacing 
accretion of hydrogen-rich material by accretion of he¬ 
lium in order to bypass the hydrogen flashes, and use 
as accretion rates the effective accretion rates obtained 
in the long-term calculations of recurrent outbursts de¬ 
scribed in (j2] These represent the actual rates of he¬ 
lium accumulation. We only consider the fraction of 
the (Mwd.AT) parameter space that leads to a Mch 
WD within a reasonable amount of time for a realis¬ 
tic donor mass. Refining the WD mass scale, we have 
performed several series of helium accretion simula¬ 
tions corresponding to three hydrogen accretion rates. 
We calculated between 500 and 1500 full consecutive 
cycles of helium accretion, followed by TNR and ejec¬ 
tion, depending on the case. 

The first evolution series is for WDs with masses 
of 1.0, 1.1, 1.25 and 1 . 34 M 0 accreting material com¬ 
posed of 98% helium and 2% heavy elements at rates 
corresponding to a constant accretion of hydrogen at a 
rate of 2 xlO“^M 0 ?/r“^. 

We found that — similarly to the hydrogen-flash 
cycles — the accreted mass per helium-flash cy¬ 
cle is higher for lower WD masses, ranging from 
^ 7 . 2 xlO“®M 0 for the I.OM 0 case to ^ 4 . 1 xlO“®M 0 
at the end of the 1 . 34 M 0 simulation. The accreted 
mass required to trigger a flash decreases gradually 
as the simulations progress; over the first couple of 
hundred cycles the decreasing rate is rapid, while later 
on it slows down. The decrease in accreted mass per 
cycle over time, at a constant rate of accretion, means 
the cycle duration decreases as well. In Fig|2]we plot 
the helium-flash cycle duration D of the growing WD. 

The ejected mass (mej) is always less than the ac¬ 
creted mass (mace) and, like the accreted mass, de¬ 
creases as the simulations progress. However, the ratio 
of ejected to accreted mass does not remain constant, 
but decreases continually and eventually reaches zero 
after a few hundred cycles, that is, there is no mass 
loss at all. We continued the simulations for several 
hundred cycles more to ascertain that this behavior — 
zero mass loss per helium flash — continues. It does. 

At the beginning of each evolution run (for a given 
initial WD mass), helium accumulates, building up a 
relatively thick helium-rich layer. Helium burns into 
heavier elements at the bottom of this layer. At first, 
the burning front advances outward in mass at a slower 
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Fig. 7.— Helium flash cycle duration (D) vs. WD 
mass (Mwd)- Different curves correspond to different 
initial WD masses; 1.0,1.1,1.25 and 1.34Mq. The 
black dotted eye-fit line demonstrates the convergence 
to an asymptotic behavior pattern. The effective accre¬ 
tion rates used correspond to accretion of hydrogen- 
rich matter at a rate of 2xlO“^MQ?/r“^. The cyan 
star represents the helium flash of cycle #2574. 


pace than helium accumulates. In time, however, over 
a few hundred cycles, the internal temperature rises, 
both at the bottom of the helium shell and through¬ 
out it. As a result, the helium burning front advances 
outwards at a higher pace and at the same time, the 
helium mass fraction throughout the helium-rich layer 
decreases, until, eventually, helium fuses completely 
into carbon, oxygen and heavier elements. Thereafter, 
helium burns at each cycle as it is accreted. 

As an example of the evolution of the WD interior 
temperature, we plot in Fig0 the temperature profile 
for the WD core of the 1.1 Mq model at a few points 
in time, marked by the cycle number and mass at that 
point. During the evolution, the temperature gradient 
throughout the star monotonically decreases. At cycle 
^5 almost the entire star is still at the initial tempera¬ 
ture, except for the outer layers which are affected (due 
to diffusion, convection, and heat conduction) by the 
intensely burning accreted helium in the envelope. As 
the evolution progresses, heat penetrates inwards and 
slowly raises the temperature deeper in the star. This is 
clearly seen in the profiles up to cycle #58; at the cen¬ 
ter of the star, however, the temperature still remains 
close to the initial value. From cycle #98 on, we see 
the center of the star heating as well. This cycle repre¬ 
sents the point in time where the helium layer becomes 



^cyc#5 (1.105M 
_^cyc#40 (1.114M 
—eye #58 (1.125M 
—eye #98 (1.157M 
-*-eye #200 (1.204M 
— eye #400 (1.247M 

0 1.0x10 “K ’X. 
ra i V. “at. 

^ St,. 

5.0x10 ^ 

1.0x10 



Fig. 8.— Profiles of temperature (top), density (mid¬ 
dle) and the Fermi parameter (bottom) at six points 
in time throughout the evolution of the I.IM© WD, 
accreting helium at rates corresponding to the con¬ 
stant rate of accretion of hydrogen-rich matter of 
2xDuring the ^6x10® years and 400 
helium flashes shown, the WD grows in mass from 1.1 
to 1.247M0, becoming much hotter and less degen¬ 
erate. The decreased degeneracy makes later helium 
flashes less violent, allowing the WD mass to grow to¬ 
wards Mch- 


depleted enough to allow quasi-steady helium burning. 
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and where the WD stops ejecting mass. 

The helium profile at the same six cycles is shown 
in Fig|9] demonstrating that the burning front is mov¬ 
ing towards the surface. The fraction of the helium in 
the accreted layers decreases monotonically, dropping 
from ^90% to ~30% over 400 helium flashes. After 
helium flash ^98 we find that helium is burnt at the 
same rate at which it is accreted. Successive flashes 
are mild and non-ejective because the surface layers 
have become hot and less degenerate. 
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Fig. 9.— Helium profiles at six points in time as de¬ 
scribed in FiglU 


The number of helium flash cycles required to reach 
the non-ejecting stage depends on the initial WD core 
temperature: the lower the initial temperature, the 
more flashes and time are needed to reach the non¬ 
ejecting stage. When the internal temperature is low, 
heat penetrates inward more easily and as a result, the 
temperature in the outer helium-rich region is too low 
to enable the accumulated helium to burn. Eventually, 
the internal temperature rises sufficiently to maintain 
the temperature in the helium-rich region above the 
helium burning threshold, thus helium burns continu¬ 
ously and the advance to the non-ejecting flashes stage 
is more rapid. This may require more than a thousand 
helium flashes. For example, the I.OMq evolution run 
began with a low initial core temperature of lO^K and 
was continued for ~800 flash cycles towards the non¬ 
ejecting stage. 

At the end of the series simulations, the WDs have 
grown significantly in mass. The I.OMq WD in¬ 
creased to 1.18Mq, the I.IMq WD became 1.26Mq, 
the 1.25Mq WD grew to 1.342Mq and the 1.34Mq 


WD became 1.382Mq. The central temperature in¬ 
creased significantly, the highest temperature, 5.5xlO*K, 
being obtained for the most massive WD (1.382Mq). 
Even more important, the increased mass in the outer 
portion of each WD is composed of carbon and oxy¬ 
gen, not helium. There is also no doubt that each of 
these WDs will continue growing in mass, and reach 
Mch, if we continue accretion at the same rate. 

In conjunction with the temperature rise in the in¬ 
terior, the WDs become more condensed (due to their 
increase in mass). Remarkably, despite the rise in WD 
densities, the WD electron degeneracies decrease due 
to rising temperatures powered by the burning helium. 
This is the central result of this paper: Carbon-oxygen 
WD masses grow as helium burning becomes steady, 
instead of violently and with ejection. Profiles for den¬ 
sity (p) and the Eermi parameter (ep = Me/kT, where 
p,e is the electron chemical potential) are presented in 
FiglDfor the same points in time as for the temper¬ 
ature profiles, showing that throughout the evolution, 
the central density increases fourfold and the Eermi 
parameter decreases threefold compared to their initial 
values. 


The second evolution series consists of WDs with 
masses of 1.2, 1.25 and 1.3Mq accreting at rates cor¬ 
responding to a constant hydrogen accretion rate of 
5x 10“^MQyr“^. As in the first series, we found that 
the accreted mass per helium cycle is higher for lower 
WD masses (ranging from ^8.0 x10“^Mq at the be¬ 
ginning of the 1.2Mq case to ~3.0x10“®Mq at the 
end of the 1.3M0 simulation), the accreted and ejected 
masses decrease gradually and the ejected mass even¬ 
tually becomes zero. At the end of the simulations 
the 1.2Mq CO WD grew to 1.356Mo, the 1.25Mo 
grew to 1.31Mq and the 1.3Mq increased in mass to 


1.36Mo. 

Eor both series, the peak WD luminosity attained 
during a helium flash is about IO^Lq and the peak 
helium shell burning temperature is about 6xlO®K, 
with little variation among the cases considered. These 
values a re in agreement with those report ed by other 
authors ( Kato et aP 19891: Jose et al. 1993 ) for helium 
burning shells. We emphasize, however, that the simu¬ 
lations presented here are the first to show the dramatic 
changes in helium burning, due to the change in struc¬ 
ture driven by heating of the underlying WD, when 
hundreds of successive flashes are self-consistently 
simulated. 


All models accreting helium, at effective rates cor¬ 
responding to the same hydrogen accretion rate, con- 


8 















Mv\/D[I^Sunl 


Fig. 10.— Same as Figl2]for effective helium accre¬ 
tion rates corresponding to accretion of hydrogen-rich 
matter at a rate of 5 xlO“^M 0 yr“^ and initial WD 
masses of: 1.2,1.25 and 1.3Mq. 


verge asymptotically to a common behavior pattern as 
shown in Figl7]and[T0l This may be taken to represent 
the continuous evolution of the WD of lowest mass up 
to the Mch limit. These figures demonstrate that ini¬ 
tial conditions do not affect the inevitable long-term 
outcome, that is, collapse. A continuous evolutionary 
simulation throughout the helium flash phase, starting 
with a I.OM 0 WD and ending at Mch would require 
not only a huge amount of time, but also a different nu¬ 
merical scheme, as the mass scale characteristic of the 
active region changes by several orders of magnitude, 
as the WD grows in mass. 


A third evolution series we simulated corresponds 
to a constant accretion of hydrogen at a rate of 
lO-^M 0 yr-i. The I.OM 0 WD accreted ~3x IQ-^Mq 
of helium each cycle and then ejected all of it. How¬ 
ever, in cases with initial WD masses of 1.1, 1.2 and 
1.32M0, masses of 10 “^— 1 O“^M 0 were accreted 
during the first cycle, and then the WDs reached tem¬ 
peratures of above lO^K, the nuclear and neutrino 
luminosities rose by many orders of magnitude and 
heavy element production began. These may corre¬ 
spond to the detonati on of a faint thermonucle ar super¬ 
nova as described by Bildsten & ShenI ( 2007h . Those 
authors calculate the relation between the thermonu¬ 
clear and dynamical timescales and deduce that for 
low accretion rates, masses higher than saO. 9 M 0 will 
accrete 2 x 10 “^— 1 O“^M 0 and undergo what they call 
a ”.SNIa”, to express that the intensity of the explosion 
is about one tenth of that of a typical SNIa. 


3.3. Accretion efficiency of the helium hashes - 
The parameter space revisited 

The key conclusion of the simulations just de¬ 
scribed is that the first few helium flashes do not rep¬ 
resent the long-term evolution of a rapidly accreting 
WD. The efficient mass-ejection behavior exhibited 
during the early phases of helium accretion cannot be 
extrapolated to later times. Doing so has led to the 
entirely erroneous conclusion that WDs cannot grow 
in mass to reach Mch via the SD channel. On the 
contrary, we have just shown that a WD can certainly 
grow to Mch while rapidly accreting hydrogen. How¬ 
ever, we must correct the extent of the parameter space 
that allows for this outcome, considering the additional 
loss of mass due to helium flashes. 

Since we have, in effect, evolved a WD virtually 
continuously from I.OM 0 all the way to 1.38M0, we 
can now estimate how much mass will be lost through¬ 
out the evolution due to the helium flashes alone, and 
then incorporate this estimate in our calculations of r 
and Ms to improve the accuracy of the limits on the pa¬ 
rameter space that can produce SNIa progenitors. We 
note that the following estimations will apply only to 
the hydrogen accretion rates that we have sampled for 
the helium accretion. 

We first attempt to estimate the amount of mass 
AM that the donor is required to transfer to the WD in 
order for the latter to grow from the initial WD mass 
to 1 . 4 M 0 . The required mass is thus defined as: 



where rj = {mace — mcj)lmacc is the mass reten¬ 
tion efficiency. For the first series (i.e., correspond¬ 
ing to an accretion rate of 2 xlO“^M 0 yr“^) p in¬ 
creases from about 7% at the beginning of the sim¬ 
ulation, at I.OM 0 to 100 % at I.I 5 M 0 , and for the 
second series (i.e., corresponding to an accretion rate 
of 5 xlO“^M 0 ?/r“^) T] increases from about 17% at 
the beginning of the simulation, at I. 2 M 0 to 100 % at 
1.28M0. Since the evolutionary calculations are not 
continuous, but piecewise, the integral in Eql2]can be 
written as a sum function yielding: 


AM^Y. 


Ami 


(8) 


We divide the entire evolution into several segments 
and thus obtain for the first series ( 2 xlO“^M 0 j/r“^), 
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starting from I.OMq: 


am., 


,, ,v_, 1.1-1.0 

-^series(I) ~ 7^ 

1.12-1.11 I 1.15-1.12 


53% 


97% 


1 . 11 - 1.1 
■ -^Wo — 

1.4-1.15 , 
100 % ' 

-7/ 


(9) 


I.8M0 


and for the second series (5x10 'MQyr ^), starting 
from I. 2 M 0 : 


AAfgej.jes^jj^ 


1.237-1.2 

17 % 


1.246-1.237 

20 ^ 


1.26-1.246 

30% 


1.274-1.26 I 1.28-1.274 , 1.4-1.28 
50% 60% 100% 


O.47M0 

( 10 ) 


These masses must be added to the donor mass 
(Mg) calculated in ij2]for the corresponding accretion 
rates and initial masses. Based on Figl3 for a I.OM 0 
WD accreting at 2xl0~'^MQyr~^, this would be a 
mass of ^O.58M0, which will yield a corrected es¬ 
timated minimal donor mass of ^2.38Mq, and for a 
I. 2 M 0 WD accreting at 5xl0~'^M 0 yr~^, this would 
be a mass of ~O.23M0, which will yield a corrected 
estimated minimal donor mass of ^O.TMq. The re¬ 
quired donor mass as a function of initial WD mass 
is presented in FigHIl which shows that it decreases 
rapidly with increasing WD mass. 


In order to ensure a secular stable mass transfer 
by Roche-lobe overflow, the dono r star would need to 
be at most as massive as the WD tergma & Fedorova 
I 990 I : iKniggel I 2 OIII: iToonen et al.l 2014 ). The black 


line in FigfTTIrepresents the points where the WD and 
donor are of equal mass, meaning that the area below 
this line will experience stable accretion. For exam¬ 
ple, a secondary of I.IM 0 initial mass transferring 
matter onto a WD of equal initial mass at a rate of 
2 xlO“^M 0 j/r“^, or onto an initially more massive 
WD at a higher rate (I. 2 M 0 and 5xlO“^M 0 j/r“^, 
respectively), will allow the WD to reach, eventually, 
the Chandrasekhar mass. Dividing the donor’s mass 
by the mass transfer rate, yields a relatively short time 
scale, ranging between 5 and 2 million years. 


This constraint on the secondary mass does not ap¬ 
ply, if the system is a symbiotic binary, where the WD 
accretes matter from its red giant companion’s wind. 
In this case, the secondary star may be more massive 
than the WD, the limitation being that the initial (main 
sequence) mass of the secondary be lower than that of 
the primary. It may well be that many binaries leading 
to SNIa are symbiotics, but this conjecture requires a 
separate feasibility study. 



Fig. 11.— Lower limit for the donor mass (Mg) re¬ 
quired to grow a WD to the Chandrasekhar mass, 
where the effects of hydrogen flashes and helium 
flashes are both taken into account. The effect of he¬ 
lium flashes alone is marked by dotted lines. As in 
Figl3 the black line corresponds to MwD=Afs, which 
defines the upper limit for stable mass transfer, assum¬ 
ing Roche-lobe overflow. This limit does not apply for 
wind accretion in a symbiotic binary. 


4. Discussion and conclusions 

4.1. Testing the helium accretion procedure 

In order to test whether we have not strayed from 
the hydrogen accretion evolutionary track by introduc¬ 
ing helium accretion as a bypass, we have resumed hy¬ 
drogen accretion at the end of the helium accretion run 
of the 1.25M0 model. Starting from a 1.34M0 WD 
with a central temperature of 10®K (the point where 
the evolution run for 1.25M0 model ended), we cal¬ 
culated more than 4700 additional hydrogen nova cy¬ 
cles. We found that these behaved very much like the 
I. 34 M 0 model described in 0 The small differences 
in D are attributed to the different core temperatures: 
the I. 34 M 0 hydrogen accretion model described in @ 
began with a core temperature of 3xlO^K, while the 
evolved model reached a central temperature of 10®K. 
The cycle duration affects the amount of mass that is 
accreted and ejected but has little effect on the net ac¬ 
creted mass. 

Throughout the evolution the hydrogen nova cycles 
are virtually identical, accreting ~3.4xlO“^M0 and 
ejecting ^1.5xlO“^M0 thus retaining ^1.9xlO“^M0 
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(~56%) at the end of each cycle. The ejected mass 
consists mostly of the accreted hydrogen with ^38% 
helium and ^2.4% heavy elements. The average 
nova cycle is ~2 years, the central temperature grows 
slowly with a total increase of less than ~0.5% 
and the maximal temperature per cycle is steady at 
^1.57xlO®K throughout the evolution. This unifor¬ 
mity is maintained throughout, except for a brief inter¬ 
ruption at cycle 2574, which is, in fact, a helium flash, 
exhibiting an entirely different behavior. 

Until cycle #2574, a total net mass of ~4.87x 
has been accreted and burnt into helium. In the 
course of the flash of the 2574th cycle, the WD ejects 
^4.48 x10“^Mq, retaining only ^3.85 x10“®Mq 
(~ 8 %) of the mass it has accreted since the beginning 
of the simulation. The ejected mass during this cycle 
consists of He (^^43.75%) and the rest (^56.25%) are 
heavy elements. This ejecta composition is typical of a 
helium nova flash, and the flash resembles the helium 
flashes produced by our helium accretion simulations 
described in Cycle #2574 of this simulation is 
marked as a cyan star in Fig 17] showing that also the 
duration of this cycle (~6 years) fits in with the he¬ 
lium accretion simulations. The luminosities and the 
maximal temperature, too, are similar to the helium 
accretion cases. The maximal temperature per cycle, 
presented in Fig [121 shows that throughout the evolu¬ 
tion Tjnax never became higher than ^1.6xlO®K — 
barely hot enough for fusing helium — meaning that 
the rate of heavy element production is slow. How¬ 
ever, during cycle #2574 the maximal temperature 
rises to ^8.4xlO®K — well above the threshold for 
helium fusion and for many heavy elements as well, 
explaining the heavy element enrichment of the ejecta. 
The effective temperature (FigfTTli during this cycle 
is significantly higher as well, reaching as high as 
~2.5xlO®K which has a peak black body radiation at 
~l.lKev (^l.lnm) meaning it would be detectable 
mostly in the soft X-ray band. The bolometric, nu¬ 
clear and neutrino luminosities are all significantly 
higher during this cycle (FigfTBll. A typical hydrogen 
cycle has a maximum bolometric luminosity, LboU 
of ~ 6 xlO^L©, and a neutrino luminosity, Tneut, of 
^I.OLq. During the helium cycle, Lboi becomes 
as high as ^5.4xlO®L©, that is, almost ten times 
brighter, and Lneut rises to ~7.5xlO^T©, nearly four 
orders of magnitude higher. 

After the helium flash, the WD relaxes for ~3 years, 
after which the hydrogen flashes are resumed. The first 
few cycles are slightly irregular, but after a short period 



Time[yr] 


Fig. 12.— Effective temperature (Tea) and maximal 
temperature (Tmax) on a logarithmic scale vs. time. 


of adjustment, the evolution continues in the same typ¬ 
ical fashion as before the helium flash. This trend can 
be seen in FigfTTl and FigfTBl where the temperatures 
and luminosities slowly return to normal over the cy¬ 
cles following the helium flash. 

In conclusion, the identical behavior of the helium 
flash of cycle #2574, resulting from the self-consistent 
simulation of 2573 preceding hydrogen flashes, and of 
the helium flashes obtained by adopting helium ac¬ 
cretion, indicates that the shortcut that replaces the 
simulation of hundreds of thousands of hydrogen nova 
cycles, does not change the simulation results. 


4.2. Summary and conclusions 


The goal of this study was to And the conditions 
required for producing a SNIa progenitor by accre¬ 
tion onto a WD in a close binary system. We thus 
examined the parameter space spanned by WD mass 
and accretion rate, with the requirement that the WD 
grow in mass, despite recurrent nova outbursts, dur¬ 
ing which some of the accreted mass is ejected back 
into space. Prev ious studies dHillman et alJ 12015 


Yaron et al.l l2005h have shown that the relevant pa¬ 


rameter space is limited by accretion rates in the range 
0.3 — 6 xlO“^M 0 ?/r“^, for the entire WD mass range, 
0.65 — 1.4M©. The effect of the initial WD core tem¬ 
perature (intrinsic luminosity) is small in these cases 
and was ignored. 


We began by simulating the evolution of WDs with 
a range of masses, accreting hydrogen-rich (solar com¬ 
position) material at given constant rates within the 
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Fig. 13.— Luminosities (bolometeric (Lboi), nuclear 
(Lnuc) and neutrino (Lneut)) on a logarithmic scale vs. 
time. 


range. The evolution was followed for several hun¬ 
dreds of nova cycles, to determine the efficiency of 
mass retention for each parameter combination. At 
this first step, we neglected the effect of eventual he¬ 
lium flashes, bound to reduce the fraction of accreted 
mass that is retained by the WD. We were thus able 
to obtain a lower limit for the time required for a WD 
of any mass to reach Mch by accretion at fixed given 
rates in the relevant range. This resulted in the first re¬ 
duction of the initial parameter space, by demanding 
that the lower time limit be shorter than the Hubble 
time. Low-mass WDs (0.65Mq), accreting at rates of 
5xl0~^MQyr~^ and lower were eliminated as realis¬ 
tic SNIa candidates. 

The next step was to estimate the required donor 
mass. If the donor is a main-sequence star and mass 
transfer is via Roche-lobe overflow, then, in order to 
obtain a steady mass transfer rate, the secondary can¬ 
not be more massive than the primary. This results in 
severe constraints on the masses of either the WD or 
the secondary, when low accretion rates are involved. 
However, if the system is a symbiotic binary, where 
the WD grows by wind accretion, this difficulty is cir¬ 
cumvented. 

The truncated parameter space required, however, 
a further reduction due to helium flashes that were ex¬ 
pected to occur at longer intervals, interrupting the reg¬ 
ular series of hydrogen flashes. In fact, several stud¬ 
ies that considered helium flashes came to the conclu¬ 
sion that a WD cannot reach Mch by accretion from 


a companion under any circumstances. These studies 
calculated, however, only a small number of helium 
flashes and extrapolated the results, on the assump¬ 
tion that the flashes will stay the same indefinitely. 
The full evolution simulation of hydrogen accretion 
through hundreds of thousands of hydrogen flashes, 
with many hundred helium flashes in between, is pro¬ 
hibitively time consuming. 

The key methodology and result of this paper is 
that we proceeded to simulate the long-term evolu¬ 
tion through helium cycles by suppressing the hydro¬ 
gen flashes and adopting helium accretion at the ef¬ 
fective rate. We first tested by a single long-term cal¬ 
culation that such a procedure would not alter the re¬ 
sults. In fact, this procedure has been used before 
in many studies. Since even these calculations were 
extremely time-consuming, we considered the entire 
range of WD masses in the truncated parameter space, 
and sampled three accretion rates from the range. We 
found that for the entire range of WD masses, the he¬ 
lium flashes become gradually less violent and even¬ 
tually settle down to mild eruptions with no mass loss 
at all. The reason for this is that the WD is secularly 
heated by the helium flashes at its surface. Thus the 
temperature of the entire star becomes so high that 
the helium burning becomes less degenerate. As the 
temperature increases, the accumulated helium shell 
is slowly depleted, and Anally the helium is steadily 
burnt at the rate which it is accreted. This occurs after 
a few hundred flashes. Thus the total mass retention 
efficiency is indeed initially diminished by degenerate 
helium flashes. But after of order 100 helium flashes 
the heated WD retains all of the matter it accretes. For 
example, a 1.085Mq secondary transferring mass at a 
rate of 2xlO“^M0yr“^ to a WD with an initial mass 
of at least I.OSSMq would permit the latter to reach 
Mch in a reasonable time. 

In conclusion, accounting for both hydrogen and 
helium flashes over the entire evolution time, we And 
that there is a significant region of parameter space 
where the single-degenerate scenario for SNIa is valid. 
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